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ABSTRACT 

We present Submillimeter Array observations of high frequency SiO masers around the supergiant 
VX Sgr and the semi-regular variable star W Hya. The J = 5 — 4, v = 1 ^*SiO and v — ^^SiO masers 
of VX Sgr are shown to be highly linearly polarized with a polarization from 5 — 60%. Assuming 
the continuum emission peaks at the stellar position, the masers are found within ~ 60 mas of the star, 
corresponding to 100 AU at a distance of 1.57 kpc. The linear polarization vectors are consistent with 
a large scale magnetic field, with position and inclination angles similar to that of the dipole magnetic 
field inferred in the H2O and OH maser regions at much larger distances from the star. We thus show 
for the first time that the magnetic field structure in a circumstellar envelope can remain stable from 
a few stellar radii out to ~ 1400 AU. This provides further evidence supporting the existence of large 
scale and dynamically important magnetic fields around evolved stars. Due to a lack of parallactic angle 
coverage, the linear polarization of masers around W Hya could not be determined. For both stars we 
observed the ^^SiO and ^^SiO isotopologues and find that they have a markedly different distribution 
and that they appear to avoid each other. Additionally, emission from the SO 85 — 44 line was imaged 
for both sources. Around W Hya we find a clear offset between the red- and blue-shifted SO emission. 
This indicates that W Hya is likely host to a slow bipolar outflow or a rotating disk-like structure. 

Subject headings: Masers, Polarization, Stars: Individual (VX Sgr, W Hya), Stars: Late- Type 



1. INTRODUCTION 

Many different SiO maser transitions have been ob- 
serve d around evolved stars of different classes ([Kemballl 
l2007l and references therein). High angular resolution 
observations of 43 and 86 GHz masers have shown that 
they occur within a few stellar radii from the photo- 
sphere, in the dynamical region bet ween the pulsating 
star and the dust for r nation zone (e.g. ICotton et alll2004l: 
iGreenhiU et allll995t IChen et al.ll2006D . Si O masers are 
therefore excellent probes of the processes that drive stellar 
mass loss, and that define the outflows that give rise to po- 
tential asymmetries in the circumstellar envelope (CSE). 

SiO maser emission from evolved stars has been detected 
from the u = to u = 4, J = 1 — up to at least the 
J = 8 — 7 transitions (e.g.. Gray, Humphreys & Yates 1999; 
Pardo et al. 1998; Humphreys et al. 1997; Jewell et al. 
1987). The conditions giving rise to emission are typically 
n(H2)=10i°=^i cm-3 and T^ > 1500 K, with the dominant 
pumping mec hanism, radiative or co l hsional, a matter 
of de bate (e.g. iHerpin fc Baudrvll2000l: iLockett fc Elitzuij 
Il992| ). However, it is agreed that line overlaps between the 
SiO main isotopologues (^^SiO, ^ssiO and ^oSiO) likely 
play a role in the pumping of e.g., the very highly excited 
V — 4: J = 5 — 4 maser lin e 7000 K above ground 
state; ICernicharo eTallliggSD . SiO maser v = 1 and 2, 
J = 5 — 4 emission is commonly found toward Mira vari- 
ables for which the lower- J masers have already been de- 
tected. High-frequency SiO maser emission in Miras ap- 
pears to be strongly tied to stellar pulsation, with the ve- 
locity of the emission a function of phase, and this phase 
dependence correlated with photon luminosity. Indeed, 
the J = 7 — 6 and 8 — 7 lines are we ak or absent from 
d> = 0.4 to 0.7, unlike the low- J masers (|Grav et al.llT99l . 
iHumphrevsl ((1999) flnds that highly rotationally excited 
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SiO maser lines arise from a subset of the physical condi- 
tions leading to J = 1 — emission, and are strongest in 
dense, warm post-shock gas. 

Besides probing the dynamics and physical conditions 
in the CSEs, masers are also go od probes of the mag- 
netic held (e.g. IVlemmingsl l2007l) . Most of the informa- 
tion on the magnetic field around evolved stars comes 
from maser polarization observations. Ordered magnetic 
fields with a strength of order a mC have been detected 
in th e OH maser regions a t large distances from the star 
(e.g. ISzvmczak et al.lll998l ). Closer in, H2O masers also 
indicate the presence of a dynamically important mag- 
netic field, with typical values of a few hundred mG (e.g. 
|Vlemmin^_eLalJ l2005f). Finally, 43 and 86 GHz SiO 
maser observations reveal field strengths of several Gauss 
at only a few stellar radii fro m the star, assuming a stan- 
dard Zeeman interpretation (jKemball fc Diamondl 119971 : 
iHerpin et al.| [2006) . The higher frequency SiO masers dis- 
play a large fractional linear polarization (jShinnaga et al.l 
l2004l hereafter S04) , and this polarization is potentially a 
good probe of the circumstellar magnetic field morphology. 
However, as these masers likely exist in the regime where 
the Zeeman splitting is comparable to the rate of stimu- 
lated maser emission, the i nterpretation of their polar iza- 
tion is not straightforward (jNedoluha fc Watsonlll994l ). 

Here we present the polarization and distribution of the 
w = land2J = 5- 4 ^SSiO and u = 0, J = 5 - 4 ^SSiO 
masers around the supergiant VX Sgr and the semi-regular 
variable st ar W Hya obser ved with the Submillimeter Ar- 
ray (SMA; iHo et al.ll2004[ ). In the remainder of the paper 
we will forgo most mentions of the rotational quantum 
numbers J = 5 — 4 and only give the vibrational quantum 
number w = 0, 1 or 2. In addition to the maser observa- 
tions, we present observations of the SO 65 — A4 line that 
was detected for both stars. 

VX Sgr is a red supergiant star, lo cated at a distance 
of 1.57 ± 0.27 kpc ()Chen et al.ll2007D . It has an optical 
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Fig. 1.— SMA spectra of VX Sgr (left) and W Hya (right) for the ^^SiO v = l,J = 5-4 masers at 215.6 GHz (thick sohd hne), the ^^SiO 
v = 2, J = 5- A masers at 214.1 GHz (thick dashed hne), the ^^SiO v = 0, J = 5 — A masers at 214.4 GHz (thin soUd Unc) as well as the SO 
55 — 44 line at 215.2 GHz (thin dotted line). The vertical lines denote the stellar velocity. 



variability period of 732 days which suggests, assuming 
the Mira period-luminosity relation, it has a mass of ~ 
10 Mq. Its envelope hosts a number of maser species, 
from which the stellar v elocity is estimated to be V|sr = 

5.3 km s-i (jChapman fc Cohen 1 986,). 

W Hya, at a distance of 98 pc ()Vlemmings et al.ll2003l) . 
is typically classified as a semi-regular variable star as its 
period is observed to vary between ^ 350 — 400 days. 
However, it is also often considered a Mira variable due 
to the reg ular shape and the large ampli tude of its light 
curve (e.g. iGomez Balboa fc Lepinelll986D . The envelope 
of W Hya harbors many different maser species, and its 
OH maser spectrum indicates a stellar velocity of Visr — 
40 km s-i (ISzvmczak et al.lll998[ ). 

2. OBSERVATIONS 

The SiO masers of VX Sgr and W Hya were observed 
on UTC 2008 July 20 with the SMA in extended config- 
uration. The resulting beam size was ^ 1.5 x 1.0 arc- 
seconds. The observations were centered at 214.8 GHz 
(A = 1.4 mm) to allow simultaneous observation of both 
the V = 1 and v = 2 ^^SiO maser lines at 215.5959 and 
214.0885 GHz respectively in one 2 GHz band. This addi- 
tionally covered the frequencies of the ^^SiO v = maser 
and the SO 55 — 44 line. The frequency resolution was 
400 kHz, corresponding to 0.56 km s"^ The 2 GHz 
bandwidth, excluding the channels with line emission, was 
used to observe the dust continuum. For the relative am- 
plitude and phase calibration of VX Sgr we used NRAO530 
(~ 12.4° away; 2.45 Jy) and for W Hya we used 1337-129 
(~ 15.8° away; 2.4 Jy). For the absolute amplitude cali- 
bration we used 3C454.3 with a flux density of 22.63 Jy. 
In both cases we used the brightest v = 1, ^^SiO maser to 
self-calibrate the data. We estimate the absolute positions 
to be accurate to ~ 0.2 arcseconds, while the relative posi- 
tional accuracy of the different maser species and the dust 
should be good to better than ^ 15 mas. The data were 
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Fig. 2. — The feed position angle (P.A.) vs. the apparent flux 
density for four of the ^°SiO d = 1, J = 5 — 4 masers around VX Sgr. 
The error bars denote the uncertainties on the flux density as derived 
by the jmfit Gaussian fitting procedure and do not contain potential 
systematic effects due to time variable errors in the gain calibration. 
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initially calibrated using the MIR IDL package, continuum 
subtraction and averaging was performed in MIRIAD, and 
the imaging and phase self-calibration was performed in 
AIPS. The channel rms noise is ^ 0.05 Jy beam~^. 

After imaging and self-calibration we used the AIPS 
image-plane component fitting task jmfit to determine the 
position of the maser emission in each individual velocity 
channel. When emission is detected at > lOcr, we define 
the fitted position i n each channel as a maser spot follow- 
ing the definition of lChen et al.l ()2006l ). A maser feature is 
then a group of maser spots that occurs within a small spa- 
tial and spectral region. Very Long Baseline Interferome- 
try (VLBI) SiO maser observations at high frequency reso- 
lution indicate that there are typically several tens of such 
features with line- widths of order 1 km s~^and a typical 
size of ~ 1 AU (e.g. lChen fc Shenll2008[ ). Both the spectral 
as well as the spatial resolution of our observations are in- 
sufficient to identify individual maser features. The maser 
spots defined here thus denote the flux weighted average 
position of features that contribute flux in the individ- 
ual velocity channels. As the beam size (half-power beam 
width, HPBW) of our observations is ^ 1.5 x 1.0 arcsec- 
onds, the formal position errors on the masers depend on 
the signal-to-noise ratio (SNR) at which the maser spots 
were detected through A (5, a) = HPBW/ (2 x SNR). 

As was performed by S04, we measured the linear po- 
larization by using earth rotation polarimetry, the diurnal 
rotation of the sky over the SMA. As the SMA anten- 
nas have alt-azimuth mounts and fixed linearly-polarized 
feeds at the Nasmyth focus, the position angle (P.A.) of the 
feeds on the sky rotates during the observations according 
to Equation 1 of S04: 

PA.n=45°-a-sinS^5i^4?^- W 

cos(a) 

Here a and h are the elevation and hour angle of the 
source respectively, and 0iat is the latitude of the SMA 
(~ 19°49'27"). For VX Sgr, h ranged from -2.3 to 2.8 hr 
and the P.A. varied from —35° to 60°. For W Hya, a de- 
layed start of the observations resulted in an hour angle 
range of 0.6 to 3.4 hr and thus a smaller P.A. coverage 
from 15° to 75°. 

3. ANALYSIS AND RESULTS 

3.1. SiO Maser Spectra 

The spectra for the SiO maser transitions as well as the 
SO line are shown in Fig. [TJ Around VX Sgr, no signifi- 
cant emission was detected for the ^^SiO v = 2 transition. 
The V = 1 spectral profile has signific antly changed com - 
pared with the detection in 1986 by iJewell et al.l ([l983) 
and the peak flux density decreased from ~ 15 to ~ 6 Jy. 
The ^^SiO emission is weaker, with a peak flux den- 
sity of 4 Jy and extends to slightly broader velocities 
with respect to the star. It also shows fewer individual 
features in the spectrum. Both lines peak close to the 
Visr = 5.3 km s~^ stellar velocity. 

The SiO maser emission of W Hya is stronger than 
that of VX Sgr in all the three lines observed. Although 
the peak flux density has decreased by a factor of 5, the 
^®SiO V = 1 spectrum i s rem arkably similar to that ob- 
served by iJewell et al.l ()1987D over two decades earlier. 
It does not, however, resemble the spectrum obtained in 



1995 (|Humphrevs et al.lll997D . The shape of the ^^SiO 

V = 2 spectrum is significantly different from that of the 

V — 1 transition. Where the v = 1 line peaks around 
V^sr — 39 km s^^ with a flux density of ~ 11 Jy, the v = 2 
maser is weaker (~ 3 Jy) and peaks at ~ 44 km s~^. As 
was the case for VX Sgr, the ^^SiO maser spectrum has 
less structure than that of the ^^SiO masers, while it peaks 
close to the stellar velocity of Vlsr = 40 km s~^. At 15 Jy, 
it is stronger than the ^*SiO lines. 

3.2. Linear Polarization Analysis 

Because of the rotation of the polarized feeds on the sky 
as described in § [2l the brightness S* of a linearly polar- 
ized source changes as a function of feed position angle 
according to the equation 

5 = 5o + 5pCos[2(P.A. -x)]. (2) 

Here 5*0 is the average and Sp the polarized brightness, and 
X the electric- vector position angle (EVPA) . The fractional 
linear polarization Pi = Sp/So. An example for four of 
VX Sgr maser spots is shown in FigH) 

For a sufficient range of P.A. it is thus possible to deter- 
mine the EVPA and fractional polarization. We performed 
a least squares fit to all channels with significant (> lOcr) 
maser emission, having split the data in single 20 min 
blocks. The flux density and associated uncertainties of 
the individual observing blocks were determined using the 
AIPS task jmfit. Assuming a constant dust continuum flux 
density during the observations, we determined an initial 
gain factor for each of the time steps. The gain correc- 
tions derived from the dust were < 6%, indicating that 
the dust was not polarized at a level higher than ~ 6%. 
This, as well as the non-detection of polarization in the 
SO line, indicates that the observed SiO maser polariza- 
tion is not an artifact of our self-calibration and further 
reduction process. Still, as using the dust for gain calibra- 
tion could in principle introduce a systematic bias due to 
possible dust polarization at a level of a few percent, we 
also adopted the method described by S04 in an iterative 
minimization routine. This yielded final gain corrections 
for the individual time intervals of 2 — 9%. Still, remain- 
ing time variable uncertainties in the flux calibration will 
due to, for example, variable phase coherence add non- 
Gaussian errors to the flux density determinations. The 
influence of remaining systematic errors on the polariza- 
tion fitting can be estimated us ing the bootstrap technique 
fe.g. lEfron fc Tibshir anil [19911 ). This technique allows de- 
termination of the the fitted values and their associated 
uncertainties without relying on Gaussian statistics and 
has previously been applied to assess image fidelity in radio 
interferometric pola rimetry (|Kemball fc Martinsekl 120051 : 
iKemball et al.|[20Tol ). The bootstrap technique consists of 
the following; we treat the N data points Si for each 20- 
min observing interval as a set of independent estimates 
for the flux density. We then generate data samples of 
size N generated from the original set with replacement, 
discarding samples with fewer than four independent mea- 
surements. We generate 100,000 of such data samples for 
which the least-squares fitting routine will determine their 
best fitting 5*0, Sp and x- From the distribution of these 
sample we then determine the most probable value for each 
of the parameters as well as their most compact 68% con- 
fidence interval. 




Fig. 3. — Positions and polarization of the d = 1,J = 5 — 4 ^^SiO (left) and d = 0, J = 5 — 4 ^^SiO masers (right) around VX Sgr. The 
masers spots are plotted with respect to peak of the dust emission (aj2000 = 18''08'"04^.068 ± 0.014 and (5j2ooo = -22°13'28".416 ± 0.200). 
The uncertainty in the dust position relative to the masers is indicated by the error bars in the center. Also note the difference in velocity scale 
between the two transitions. The symbols are linearly scaled according to brightness (^*SiO; 3.8-6.9 Jy beam~^; ■^^SiO: 2.3—4.2 Jy beam~^). 
For -^^SiO, we only present the masers that were detected with a SNR > 50. For the weaker ^^SiO we plot the masers detected at SNR > 30. 
The minimum and maximum positional errors for the plotted maser spots are indicated by the error bars in the bottom left. The black vectors 
are the observed polarization vectors scaled linearly according to polarization fraction. The long dashed inner circle indicates the star with 
a radius of ~ 4.35 mas (JMonnier et al. 2004) and the solid circle indicates the location of the 43 GHz SiO masers (Chen et al. 200^. The 
dashed lines indicate the position angle and its uncertainty of the i nferred orientation of the dipole magnetic field of VX Sgr observed using 
H2O and OH masers l lVlemmings et al.|[2005l: ISzvmczak et al.ll200ll) . 



3.2.1. VX Sgr 

The polarization fractions and angles for the SiO masers 
of VX Sgr are given in the online Table. 1 and shown in 
Fig. [21 The uncertainties are derived using the aforemen- 
tioned bootstrap method. We find that the average polar- 
ization fraction of the ^^SiO v — 1 masers (Pi) = 26± 16% 
with the highest value reaching almost 80%. The maser 
EVPA X has, as seen in Fig IH a bimodal distribution, 
with peaks at (x) = -44 ±13° and (x) = 29 ± 11°. 
The ^^SiO masers have a slightly lower average polariza- 
tion {^^Pi) = 22 ± 12% with the highest fractional po- 
larization measured being 60%. The maser EVPA dis- 
tribution is somewhat wider and not bimodal, and has 
(^^x) = ~47 ± 22°. No significant polarization was mea- 
sured for the SO emission. 

3.2.2. W Hya 

The P.A. range for W Hya was too small to perform 
independent fits to both EVPA and polarization fraction. 
We thus performed a number of fits with fixed Pi or x. 
However, these did not allow for strong conclusions beyond 
that the mean fractional linear polarization is consistent 
with those found for the masers fo VX Sgr. 

3.3. Maser Distribution 

The maser distributions of VX Sgr and W Hya are 
shown in Figs. |3] and [S] The formal position errors on 
the masers are large and depend on the signal-to-noise 
ratio (SNR). We therefore only present those masers de- 
tected with a high SNR, whose positions are determined 
with better than ~ 10 mas accuracy. The distributions of 
the ^^SiO and ^^SiO masers around VX Sgr are presented 
individually with respect to the peak of the dust emis- 
sion, which was determined with an accuracy of ~ 14 mas 



(corresponding to ~ 24 AU). 

Some evidence for asymmetry in the dust envelope 
of VX Sgr has been seen using infrared interferometry 
(|Monnier et al.|[200^ . and also the dust emission peak of 
the super giant VY CMa appe ars offset from the SiO maser 
emission (jMuller et al.|[2007t ). However, if the dust peak 
of VX Sgr denote the stellar position, the ^^SiO maser 
distribution is quite asymmetric. The ^^SiO maser spots 
however, do appear to occur around the dust peak. For 
W Hya no clear offset between the dust and the masers 
was seen, and the positions for all three maser transitions 
are presented together with respect to the peak of the dust 
emission. Still, because of possible asymmetry in the dust 
emission this is not necessarily the stellar position. 

3.4. SO Emission 

Emission from the SO 65 — 44 transition at 215.2 GHz 
was detected for both VX Sgr and W Hya and is shown 
in Fig. [T] Although more typically observed at higher 
and lower frequency transitions, sulfur bearing species are 
frequently found i n the envelopes of 0-rich stars (e.g. 
lOmont et al.lll993h . Sulfur is very reactive with OH and 
will form SO. The SO emission of both VX Sgr and W Hya 
covers only a small velocity range around the stellar ve- 
locity, a range comparable to that of the SiO masers. In 
both stars the emission peaks near the stellar velocities 
with a fiux density of ~ 1 Jy. For W Hya, we find that, as 
indicated in Fig. [71 the peak of the integrated red-shifted 
emission between 40 — 47 km s~^ is offset from that of 
the blue-shifted emission between 33 — 47 km s~^. The 
separation between the emission peaks is 0.29 ± 0.04 arc- 
seconds, which at the distance of W Hya corresponds to 
~ 28 AU on the sky. The offset is almost exactly in the 
North-South direction, with a position angle of 3° ± 10. 
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Fig. 4. — The distribution of the electric- vector position angles (x) 
for the D = l, J = 5 — 4 ^*SiO masers of VX Sgr. The dashed lines 
indicate the error weighted average peaks of the bimodal distribution 
of X- The histogram was constructed using all detected maser spots 
indicated in Tabid. 
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Fig. 5. — Positions of the J = 5 — 4 transitions of ^^SiO 
V = and ■^^SiO v = 1,2 around W Hya. Positions of the 
masers are with respect to the peak of the dust emission (aj2000 = 
13''49™01^953 ± 0.015 and <5j2ooo = -28°22'06".208 ± 0.200) for 
spots with a SNR> 20. The dust position is given in the center with 
positional uncertainties. The symbols are scaled linearly according 
to brightness, ranging from 2.2-11.1 Jy beam~^. The minimum and 
maximum error bars of the plotted SiO maser spots are shown in 
the bottom right corner. The outer, dashed, circle indicates the in - 
ner dust shell as determined from models bv lWishnow et ahl 1120101 ). 
The solid lines indic ate the radius of the 43 G Hz SiO maser observa- 
tions at two epochs (|Cotton et al.l2004l . l200Sf l . The dashed ellipse is 
the el ongated radio-photosphere as determined by IReid fc MentenI 
0003). 



The SO emission of VX Sgr does not show a significant 
velocity separation and peaks at the same position as the 
dust as shown in Fig. |6l 

3.5. Dust Continuum 

The A = 1.4 mm continuum emission toward both W 
Hya (Figlll) and VX Sgr is unresolved. The integrated flux 
density of the W Hya emission is Si, = 213±3 mJy, a value 
comparable with previous measurements of 270 ± 15 mJy 
(SMA 1.1 mm; MuUcr et al. 2008), 280 ± 17 mJy (SIMBA 
1.2 mm; [Dchaes ct al. 2007), and 280 ± 30 mJy (JCMT 
1.1 mm; [van der Veen ct al. 1995), taking into account 
a lower stellar contribution to the flux at longer wave- 
length. In the case of VX Sgr, the total flux density 
Si, = 95 ±2 mJy. To our knowledge no previous millimeter 
dust measurements of VX Sgr are available for compari- 
son. 

4. DISCUSSION 

4.1. The Location of the SiO Maser Transitions 

The maser distributions shown in Figs. [3] and [5] do 
not show the clear circular morpho logy often seen in 
the 43 and 86 GHz S iO masers fe.g. ICotton et all 120041: 
ISoria-Ruiz et al.1 [20041 ) . However, with the angular reso- 
lution provided by the SMA we cannot reach the sub-mas 
positional accuracy obtained with very long baseline inter- 
ferometry for the lower frequency transitions. The ^*SiO 
maser spots of VX Sgr are offset from the peak of the 
dust emission by 28 ± 14 mas towards the North- West. 
However, the ^^SiO masers do show a hint of a circular 
distribution with an inner radius of 18 mas, slightly 
more than the ^ 12 mas observed for the 43 GHz masers 
(|Chen et al.ll2006[) . The situation for W Hya is less clear, 
although for this source 43 GHz mas ers also display onl y 
an incomplete ring morphology fe.g. ICotton et al.ll2008[ ). 
Our observations do let us conclude that specifically the 
strongest of the ^®SiO and ^^SiO masers avoid each other, 
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Fig. 6. — Contour plots of the SO Ss — 44 emission around 
VX Sgr. The contours are drawn at 10% levels of the peak bright- 
ness (795 mJy beam-^). The cross denotes the peak position of the 
dust emission. 
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something that was also found for the 43 GHz ^ ^ SiO a nd 
29SiO masers of IRC+10011 bv iSoria-Ruiz et alJ (pOOl . 

4.2. Linear Polarization 

Masers can become li nearly polarized under the influ- 
ence of a magnetic field (|Goldreich et al.lll973[ ). Here we 
measured the linear polarization of the J = 5 — 4 SiO 
masers of VX Sgr to be up to ~ 80%. This is consis- 
tent with the linear polarization fractions measured for 
the w = 1, J = 5 - 4 28SiO masers of VY CMa in S04. For 
the J = 1 — SiO maser transitions, such polarization frac- 
tions can be reached du e to interaction with a m agnetic 
field of a few Gauss (e.g. [Western fc Watsonlll98"3 ). How- 
ever, for the higher J transitions, fractional linear polar- 
ization arising solely from interaction of a saturated maser 
with a magnetic field is expected to be much lower. Al- 
ready for the J = 2 — 1 SiO maser transition, similar po- 
larizatio n fractions would need m agnetic field strengths of 
> 10 G ([Western fc Watsonlll984l ) and for J = 5-4 transi- 
tions this will be even higher. However, in the presence of 
anisotropic pumping, the linear polarization fraction can 
actually incr ease with the angular mom entum of the in- 
volved state (|Nedoluha fc WatsonlfTool) . 

The high linear polarization is therefore largely due 
to anisotropic pumping of the maser. We now need 
to assess if the EVPA can still be used to derive the 
magnetic field morphology . As described in for example 
iNedoluha fc WatsonI ()1990D this depends on the ratios be- 
tween the maser decay rate F, the maser stimulated emis- 
sion rate R and the Zeeman frequ ency shift qO. Assu ming 
a typical magnetic field of 3.5 G ()Herpin et al.ll20Q6h . the 
Zeeman frequency shift of the J = 5 — 4 SiO masers is 
^ 800 s~^. The maser decay rate F ~ 5 s"^. The rate for 
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Fig. 7. — Contours of the integrated red-shifted (Viai. > 40 km s~^; 
thin red) and blue-shifted {Vi^^ < 40 km s~^; thick blue) SO 
55 — 44 emission around W Hya over-plotted on the dust emission 
(greyscale). Contours are drawn at intervals of 0.1 Jy beam~^, start- 
ing at 0.05 Jy beam~^ (^c). The star denotes the peak position of 
the dust emission. 



stimulated emission of the J = 5 — 4 SiO masers is given by 
i? w 4 X 10-^ Tb An (S04), where Tf, is the maser bright- 
ness temperature and Ail its beaming angle. Assuming a 
maser spot size of ~ 0.4 AU similar to that of the VX Sgr 
43 GHz SiO masers (jGreenhill et al.lll995D . the brightest 
masers in our observations have ~ 2 x 10^ K and the 
weakest ones have a Tt, almost two orders of magn itude 
less. Taking Aft ~ lO'^ sr fe.g. lReid fc Moranlll988D . this 
yields R ^ 80. Therefore qV, > RS- For these values, 
the results from INedoluha fc WatsonI (jl990t ) indicate that 
the linear polarization still traces the magnetic field direc- 
tion, either parallel or perpendicular, even in the presence 
of anisotropic pumping. Specifically, the polarization vec- 
tors are parallel to the magnetic field if the angle between 
the field and the line-of-sight 9 < 9„it ~ 55° and perpen- 
dicular if > ^criti w here 6'crit is the "van Vleck angle" 
(iGoldreich et al.lll973D . 



4.3. The Shape of the Magnetic Field around VX Sgr 

Observations of 1612 MHz OH maser polarization sug- 
gest that, in the OH maser region at ^ 1400 AU from 
VX Sgr, the magnetic field has a dipole shape with a posi- 
tion angle 8 — 210°±20 and an inclin ation angle i between 
20 — 30° from the pla.ne of the sky (jSzvmczak fc CohenI 
Il997t ISzvmczak et al.|[200lD . A similar dipole configura- 
tion, with = 220° ± 10 and i = 40° ± 5 was found from 
H2O polarizat ion observations betwe en ~ 100 — 300 AU 
from the star (|Vlemmings et aLll2005fl . 

We have now been able to probe the magnetic field struc- 
ture at only ~ 20 — 100 AU from the star using SiO maser 
polarization observations. As seen in FiglSl the masers 
indicate a large scale structure in the field with a pre- 
ferred axis at ^ 30° or ^ —45°, which, within the quoted 
uncertainties, is consistent with a 90° change of EVPA 
direction. Such a 90° flip has previously also been ob- 
served in the 43 GHz SiO masers, and is attributed to a 
change of the angle between the magnetic field and th e 
line-of-sight 9 through 6'crit (jKembaU fc Diamond Il997t ). 
This is consistent with a magnetic field with an inclina- 
tion axis close to ^crit, as the curvature of the magnetic 
field near the star would, close to the plane of the sky 
where most SiO masers are tangentially amplified, pro- 
duce masers with 9 both somewhat larger and smaller than 
flcrit- Our SiO maser observations therefore suggest that 
the magnetic field close to the star has a morphology with 
an inclination axis i ~ 90 — ^crit = 35° and a position 
angle 9 - 135° ± 20 or 9 ~ 210° ± 20. These values are 
consistent with the likely dipole shaped field further out in 
the envelope, although a full 3-dimensional reconstruction 
of the field in the SiO maser region is complicated by the 
unknown location of the masers along the line of sight. 
Still, the SiO masers offer further strong evidence for a 
large-scale, most probably dipole-shaped, magnetic field 
around VX Sgr and, considering the other maser observa- 
tions, a position angle of ~ 210° is most likely. Combining 
information from all three maser species at 20 — 1400 AU 
distance from VX Sgr, we conclude that VX Sgr likely has 
a dipole magnetic field with 8 = 217° ± 7 and i = 37° ± 9. 

Having confirmed the existence of a similar large scale 
magnetic field from the SiO maser region out to the OH 
maser regions, we can extrapolate the fields measured on 
the OH and H2 O masers back t o the SiO maser region at ~ 
50 AU (jVlemmings et al.ll2005[) . As, for a dipole magnetic 
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field, B oc R^^ and as the OH and H2O maser distances 
from the star are only roughly known, the uncertainty on 
this extrapolation is large. Still, we find the field in the 
SiO maser region to be in the range of B ^ 25 — 100 G. 
While larger then the typical field found bv iHerpin et al.l 
(|200Ql , this is comparable to th e -B = 87 G field measu red 
in the SiO region of VX Sgr bv lBarvainis et alj (|1987| ). 

4.4. The SO Envelope of W Hya 

The detection of spatially offset red- and blue-shifted 
SO emission around W Hya reveals an asymmetry 
in the CSE. Previous high angular resolution SO ob- 
servations indicate that the SO emission can be en- 
hanc ed in coUimated outflows and/or equatorial structures 
(e.g. iDinh-V.-Trung et al] 120091: ISanchez Contreras et all 
l2000f l. This implies that W Hya harbors either a slow bipo- 
lar outflow or a rotating and possibly expanding equatorial 
disk. The position angle of the SO structure, measured be- 
tween the peaks of the red- and blue-shifted SO emission, 
is 3° ± 10, nearly perpendicula r to the strongly ellipt ical 
radio-photophere measured by iReid fc Menten (|2007f ) to 
have a position angle of 83° ± 18 and a size of 69 x 46 mas. 
Additionally, the red- and blue-shifted separation in SO 
is similar, both in scale, velocity and in position angle, 
to th at observed in the OH maser region f Szvmcza k et al.l 
Il998t ). From this it was inferred that the OH masers orig- 
inated in a weak bipolar outflow or a tilted circumstellar 
disk. Our observations of the SO emission thus further 
confirm this hypothesis, although we also cannot discrim- 
inate between an outflow, a disk or other more complex 
kinematic structure. This will have to wait for higher 
angular resolution observations with e.g. ALMA. 

4.5. Circumstellar Dust 

4.5.1. W Hya 

From the observed, unresolved, dust continuum of 
W Hya, we can determine whether stellar black body 
or thermal dust emission dominates the millimeter con- 
tinuum. We estimate the stellar black body contribu- 
tion using t he following parameters for W Hya: distance, 
I? = 98 pc (IVlemmings et al.l 20031) : stellar temperature. 
Toff = 2500 K (jJusttanont et al.ll2005h : and from this stel- 
lar radius, i?star = 3.4 x 10"'^'^ cm. Noting that Toff and 
i?star vary signiflcantly during the stellar p ulsation cycle 
and h ave uncertainties of 25 — 50% (e.g. iDehaes et al.l 
|2007| ). so that this calculation is quite approximate, we 
estimate 5'sta,r.i 4 r nm = 158 mJy. Thus, as was also found 



bv iMuller et aH ^QM) at 1.1 mm, a significant compo- 
nent of the millimetre continuum is due to the star itself 
rather than dust. Using the Rayleigh- Jeans approxima- 
tion to calculate d ust mass, M^ , assu ming optically thin 
emission, following iMuller et aTl (|2007[ ) 



2c^D'^adPd{Stot - Sstar) 



(3) 



where grain size Od = 0.2/im, grain mass density pd = 
3.5 g cm-3 and emissivity = 5.65x 10-''(i//274.6 GHz), 
typical assumed values for oxygen-rich stars. Adopting a 
beam-averaged dust temperature Td — 700 K, we obtain a 
dust mass Md — 4.9 x 10~^ Mq within a radius of about 
70 AU. 



The millimeter dust observations do not probe aver- 
aged mass loss history over the entire circumstellar en- 
velope, but instead can give a good indication of re- 
cent mass loss rate. For a W Hya exp ansion velocity 
of T4xp = 7 km s~^ (|Dehaes et al.l |2Q07| ). and a beam 
of ~ 1.5 arcseconds, we are probing material that has 
been ejected only over the past < 50 years. We derive 
a gas-|-dust mass loss rate of > 9.8 x 10~^ MQyr~^, as- 
suming a gas to dust ratio of 100. This is comparable to 
the distance-ad justed rate of 1.7 x 10~^ M0yr~^ from ISO 
observations of IZubko fc: Elitziiil (|2000[ ). and only slightly 
higher than the distance adjusted 4 x 10 "'' Mfr^yr"-^ value 
determined from ODIN observations by lJusttanont et al.l 

4.5.2. VX Sgr 

In order to estimate the central star black body contri- 
bution to the cont inuum of VX Sgr, we adopt D = 1.57 kpc 
Chen et al.l2007l). J? „t.. = LOxlO^cmandTes = 3200 K 
Monnier et al.l 2004^ . At A = 1.4 mm, we estimate that 
the stellar component is only S'star,i.4mm = 6.7 mJy. Us- 
ing Eqn. O the mass of the dust enclosed in a region of 
radius - 1300 AU is therefore Md = 2.1 x 10"" Mq. 
Here we have also assumed a beam-averaged Td — 700 K, 
based on results of[p anchi ct al. (1994). For the expansion 
velocity, iMurakawa et al.l (j2003ll find velocities of 10 and 
20 km s~"'^at the inner and outer edges of the water maser 
zone at 100 and 300 AU respectively. Adopting an aver- 
age Vcxp — 15 km s""'^, the observed region corresponds to 
an age of about 400 years. We estimate a total gas-|-dust 
mass of 2.1 x 10~^ Mq in this region, and therefore a 
recent mass-loss rate of 5.3 x 10~^ MQyi~^. This is in 
reasonable agreement wi th other mass loss measurements 
e.g. 3.2 x 10^5 Mgyr-i (|Netzeilll989D . 

5. CONCLUSIONS 

We have determined the magnetic field morphology at 
a few stellar radii from VX Sgr using SiO masers, and 
find that it is consistent with the interpretation in terms 
of a dipole magnetic field from OH and H2O maser ob- 
servations. This is the first evidence that circumstellar 
magnetic field morphology is conserved from close to the 
star to the outer edge of the CSE. As the magnetic field 
strengths measured in the OH and H2O maser regions im- 
ply the magnetic field is dynamically important, a further 
extrapolation of this dipole field to the star (cx R~'^) in- 
dicates that also in the SiO maser region the magnetic 
field dominates the kinetic and thermal energies, with an 
estimated field strength of S ^ 25 — 100 G. 

Our images of the SiO masers around VX Sgr and 
W Hya show that the two v = 1,2 rotational ^^SiO tran- 
sitions as well as the ^^SiO v — transition avoid each 
other. While a ring structure is seen in the ^^SiO line, the 
masers do not seem to be confined to a narrow region. 

Finally, our SO 65 — 44 observations reveal the possible 
presence of a slow bipolar outflow or a rotating, disk-like 
structure around W Hya. This structure is nearly per- 
pendicular to the elliptically extended emission previously 
detected from the radio-photophere. 

High angular resolution observations of submillime- 
ter SiO masers and their polarization will become rela- 
tively straightforward with the Atacama Large Millime- 
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ter/submillimctcr Array (ALMA). Extremely interesting 
science targets in their own right, the strength and high 
polarization fraction of the SiO masers also make them 
good ALMA polarization calibrators for a variety of fre- 
quency bands. 

This research was supported by the Deutsche Forschungs- 
gemeinshaft (DFG) through the Emmy Noether Research 
grant VL 61/3-1. 
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